Conformal Dynamics of Fractal Growth Patterns Without Randomness 
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Many models of fractal growth patterns (like Diffusion 
Limited Aggregation and Dielectric Breakdown Models) com- 
bine complex geometry with randomness; this double diffi- 
culty is a stumbling block to their elucidation. In this pa- 
per we introduce a wide class of fractal growth models with 
highly complex geometry but without any randomness in their 
growth rules. The models are defined in terms of determinis- 
tic itineraries of iterated conformal maps, automatically gen- 
erating the function <f>( n '(ui) which maps the exterior of the 
unit circle to the exterior of an n-particle growing aggregate. 
The complexity of the evolving interfaces is fully contained 
in the deterministic dynamics of the conformal map 
We focus attention to a class of growth models in which the 
itinerary is quasiperiodic. Such itineraries can be approached 
via a series of rational approximants. The analytic power 
gained is used to introduce a scaling theory of the fractal 
growth patterns. We explain the mechanism for the fractality 
of the clusters and identify the exponent that determines the 
fractal dimension. 



I. INTRODUCTION 

In this paper we introduce a new class of fractal growth 
patterns in two dimensions, constructed in terms of the 
conformal maps from the exterior of the unit circle to the 
exterior of the growing cluster. Until now most of the in- 
teresting fractal growth models included randomness as 
an essential aspect of the growth algorithms. Foremost 
in such models has been the diffusion limited aggregation 
(DLA) model that was introduced in 1981 by T. Witten 
and L. Sander Q. This model has been shown to under- 
lie many pattern forming processes including dielectric 
breakdown [2j, two-fluid flow ||, and electro-chemical 
deposition B ] . The algorithm begins with fixing one par- 
ticle at the center of coordinates in d-dimensions, and 
follows the creation of a cluster by releasing random walk- 
ers from infinity, allowing them to walk around until they 
hit any particle belonging to the cluster. Upon hitting 
they are attached to the growing cluster. The growth 
probability for a random walker to hit the interface is 
known as the "harmonic measure" , being the solution 
of the harmonic (Laplace) equation with the appropriate 
boundary conditions. The DLA model was generalized 
to a family of models known collectively as Dielectric 
Breakdown Models, in which the density of growth prob- 



ability is the density of the harmonic measure raised to 
a power 77 [gj . For rj = 1 one regains the DLA model; the 
interval rj £ (0, 00) generates a family of growth patterns 
from compact to a single needle. For rj = one obtains 
a growth probability that is uniform for all boundary 
points. This is known as the Eden model that was in- 
troduced originally to describe the growth of cancer cells 

The fundamental difficulty of all these models is that 
their mathematical description calls for solving equations 
with boundary conditions on a complex, evolving inter- 
face. It is therefore advantageous to swap for a simple 
boundary, like the unit circle, and to delegate the com- 
plexity to the dynamics of the conformal map from the 
exterior of the unit circle to the exterior of the growing 
cluster. For continuous time processes this method had 
been around for decades and had been used exten- 
sively. For discrete particle growth such a language was 
developed recently [p|-[TTl, showing that DLA in two di- 
mensions can be grown by iterating stochastic conformal 
maps. In this paper we employ this new language to 
define models in which the stochasticity is eliminated al- 
together, to create deterministic iterations of conformal 
maps with very interesting fractal growth properties. It 
is stressed below that these new models and their inter- 
esting properties are natural extensions of the discrete 
conformal dynamics; it may be very difficult to study 
such models with the traditional techniques in physical 
space. 

A central thesis of this paper is that the growth models 
introduced below are simpler to understand than DLA, 
even though the fractal geometry exhibited does not seem 
simpler. Indeed, we present below some tools and con- 
cepts that allow us to explain why the growing cluster is 
fractal. We present a scaling theory of the growing clus- 
ters, and identify the exponent that determines the frac- 
tal dimension. In sect. 2 we review the basic ideas of con- 
formal dynamics as a method to grow DLA and related 
growth patterns. In Sec. 3 we make the point that within 
this framework randomness can be eliminated from the 
discussion without changing the properties of the frac- 
tal growth: one can have deterministic growth rules with 
clusters that are indistinguishable from DLA. In Sec. 4 
we introduce fractal growth patterns that are obtained 
from quasiperiodic itineraries of iterated conformal maps. 
These itineraries are characterized by a winding num- 
ber W. The growing clusters have complex geometries 
and a difference appearance for every W. We propose 
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nevertheless that all the quadratic irrationals belong to 
the same universality class, and that the dimensions of 
their clusters are the same. In Sec. 5 we consider rational 
approximants P/Q to the quadratic irrational winding 
numbers W's. With rational approximants the growth 
patterns crossover from a fractal phase of growth to a 1- 
dimensional star-like growth pattern. We argue that the 
analysis of the crossover as a function of Q provides us 
with a scaling theory, allowing the introduction of uni- 
versality classes and the achievement of data collapse. 
In Sec. 6 we elucidate the mechanism for crossover from 
fractal to 1-dimensional growth, and identify the expo- 
nent that determines the fractal dimension. In Sec. 7 we 
summarize and offer final remarks regarding the avail- 
ability of a renormalization group treatment and of the 
road ahead. 



II. DISCRETE CONFORMAL DYNAMICS FOR 
FRACTAL GROWTH PATTERNS 

The basic idea is to follow the evolution of the confor- 
mal mapping <f>( n \w) which maps the exterior of the unit 
circle e in the mathematical w-plane onto the comple- 
ment of the (simply-connected) cluster of n particles in 
the physical z-plane ]9|-|Tll|. The unit circle is mapped 
to the boundary of the cluster which is parametrized by 
the arc length s, z(s) = $( ,l) (e ie ). This map <E> (n) (u;) is 



made from compositions of elementary maps 



$W( w )=$^-^ AnA >)) 



(1) 



where the elementary map (p\,e transforms the unit cir- 
cle to a circle with a "bump" of linear size vA around 
the point w — e 10 . Accordingly the map & n \w) adds 
on a new bump to the image of the unit circle under 
The bumps in the z-plane simulate the ac- 
creted particles in the physical space formulation of the 
growth process. The main idea in this construction is to 
choose the positions of the bumps 9 n and their sizes y/X^ 
such as to achieve accretion of fixed linear size bumps 
on the boundary of the growing cluster according to the 
growth rules appropriate for the particular growth model 
that we discuss. 

As an example consider DLA. In z-space we want to ac- 
crete particles according to the harmonic measure. This 
means that the probability for the nth particle to hit 
a boundary element ds equals P(s)ds, where P(s) (the 
density of the harmonic measure P,0,n2|) and ds are: 



P(s) 



1 



|$(n-l) ( e ifl)| 



(2) 
(3) 



Here e is the preimage of z(s). Accordingly the proba- 
bility to grow on an interval dQ is uniform (independent 
of 9). Thus to grow a DLA we have to choose random 
positions 9 n , and A„ in Eq.([|) according to 



A, 



Ao 



(4) 



This way we accrete fixed size bumps in the physical 
plane according to the harmonic measure. The elemen- 
tary map 0a,s is chosen as B 



(1 + A) 
2w 



1 + — 



(l + w) 



2 1 - A 



1/2' 



4>x,e(w) = e (f>\ t o(e~ i6 w) , 



(5) 
(0) 



The parameter a is confined in the range < a < 1, 
determining the shape of the bump. We employ a = 
2/3 which is consistent with semicircular bumps. The 
recursive dynamics can be represented as iterations of 
the map 0A„,e„(w), 



(7) 



The DLA cluster is fully determined by the stochastic 
itinerary {9i}" =1 . In Fig. [I] we present a typical DLA 
cluster grown by this method to size n = 10 5 . The main 
point of this paper is that the same method can be now 
used to grow a large variety of interesting fractal shapes, 
but without any randomness in the growth algorithm. 



III. DLA-LIKE CLUSTERS WITHOUT 
RANDOMNESS 

As a first example of a new model we will remove the 
stochasticity of DLA, leaving the growth characteristics 
unchanged. To this aim consider an itinerary 

9 n+1 = 29 n mod27r , (8) 

together with Eq.(^). Such an itinerary, although deter- 
ministic, is chaotic (in fact Bernoulli, Kolmogorov and er- 
godic), covering the unit circle uniformly, with <5-function 
correlation between consecutive 9 values. Accordingly, 
we expect the growing cluster to be indistinguishable 
from a DLA, as is indeed the case, cf. Fig. |2|. 

One advantage of the present formalism is that such 
a statement can be made quantitatively, not by eyeball. 
The function $( n '(w) and <p\^(w) can be expanded in a 
Laurent series in which the highest power is w |9|JlCt|: 



W - 



F, 



{n) 



I<_ 1 W 



(9) 



The recursion equations for the Laurent coefficients of 
can be obtained analytically, and in particular 
one shows that MuM 



in) _ 



U^ 1 + Xk 



(10) 



k=l 
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The importance of this lies in the fact that deter- 
mines that fractal dimension of the cluster. Defining R n 
as the minimal radius of all circles in z that contain the 
n-cluster, one can prove that fll3|] 



Rn < 4Fi 

Accordingly one expects that 



(n) 



(n) 



Ao , 



(11) 



(12) 



as y/Ao is the only length scale in the problem. We can 
thus present, as an example, plots of F± for our deter- 
"|) together with in any stochastic 
Another comparison is furnished 
For the DLA case it was shown 



ministic model 
DLA growth, see Fig 
by the statistics of A n 
in @0 that 

(A„ 



aDn 



(13) 



where the average is taken over the harmonic measure. 
This is in agreement with the "electrostatic relation" de- 
rived by Halsey, [ {T^j . In the Bernoulli itinerary there is 
no randomness and no probability measure, but we may 
still define a "running average" by, say, the last M iter- 
ations 



(A* 



I M 



1 

M 



E 



(14) 
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In Fig. |we show a related quantity, (X^ln-M k\k)/M 
for M = 1000 and M = 10000. We see that up to the 
expected fluctuations it settles down very quickly on the 
appropriate value of the DLA cluster, i.e. 1/aD — .877... 
Any other quantitative comparison that one can think of 
leads to the same conclusion, i.e the Bernoulli itinerary 
is a bona fide generic DLA. Of course, this is not sur- 
prising: the correlation properties of successive values of 
8 n in (||) are indistinguishable from random numbers on 
the interval [0, 2tt]. Nevertheless, our point is that the 
present growth algorithm gives us freedom to choose de- 
terministic itineraries resulting in DLA or other growth 
patterns, and we next exploit this freedom to explore new 
geometries. 



If we choose W rational, W = P/Q, then after a cross- 
over time the cluster grown is locked into a 1-dimcnsional 
object made of rays. In the next subsection we present 
an extensive discussion of the cross over time and of the 
properties of the 1-dimensional phase of growth. As an 
example consider in Fig. [5] the cluster resulting from jl5| ) 
with W — 233/144. On the other hand, for an irrational 
winding W the itinerary is ergodic and the cluster grown 
is geometrically non-trivial. As a first example we present 
the case W = p where p is the Golden Mean, p = (y/5 + 
l)/2. The fractal cluster that is associated with this rule 
is shown in Fig.|[ The cluster has a fractal dimension 
D = 1.86 ± 0.03, as determined from the scaling of . 
This is considerably higher than DLA (for which D w 
1.71). 

The golden mean is best approximated by the contin- 
ued fraction representation 



(16) 
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Such a continued fraction is denoted below as [0,1]. 
It is known that the golden mean is special in 
presenting the slowest converging continued fraction. 
Other quadratic irrationals also have periodic contin- 
ued fractions that converge faster. In Figs. [7|-|ll| we 
show the clusters grown with W — \/2, y/3, (1 + 
\ZI6)/3, (a/13 - l)/2 and \/7 respectively The contin- 
ued fraction representations of these winding numbers are 
[1,2], [1, 



1,~2], [1,2,1], [1,3], [2, 1,1, 1,4] respectively. In 
choosing these examples we picked quadratic irrationals 
whose representations converge relatively slowly. This fa- 
cilitates the exposition of scaling theory presented below. 

We note that the clusters shown have very compli- 
cated geometry. Consider for example the cases W — 
(\/l3 - l)/2 and W — \fl shown in Figs. [h| [n] re- 
spectively. They exhibit thin spiral growth patterns at 
their root, and then become bushy and thin in an ap- 
parently oscillatory fashion. Accordingly, it becomes un- 
clear whether the different quadratic irrational winding 
numbers result in the same overall fractal dimension. 
This question warrants some extra analysis. We will ar- 
gue below that in spite of the difference appearance and 
the oscillations in the "bushiness", the clusters grown 
by quadratic irrational winding numbers have the same 
fractal dimension D. 



IV. FRACTAL GROWTH WITH 
QUASI-PERIODIC ITINERARIES 

A. Winding numbers and geometry 

A new class of models is obtained by using a quasi- 
periodic itinerary. Consider a simple map of the circle 
with a winding number W: 

i+ i = 0, + 2ttW . (15) 



B. Different growth rules: period doubling itinerary 

Clearly, one can come up with an arbitrary number 
of different growth rules. In this paper we will consider 
only one additional itinerary, to underline the fact that 
quadratic irrational windings lead to a class of their own. 
This itinerary is constructed such that after every 2™ it- 
erations the points 9k chosen on the circle are equidis- 
tributed without repetitions. The order of visitation is 
determined by the following rule: 
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Qi = 2-KXi 



x l+ i = Xi + 



ki = -[log 2 (1 - x 



- 1 



(17) 



where [. . .] stands for the integer value. We refer to this 
itinerary below as the "period doubling" algorithm. The 
cluster grown with this rule is shown in Fig. ^ . The 
dimension of this cluster is D — 1.77 ± 0.02. In contrast 
to the quadratic irrationals in this case a comparison of 



F± of this cluster to F^ n> of the golden mean itinerary 
shows a different scaling dependence on n (cf. Fig. [l^ ). 



An) 



C. Universality classes? 

In the previous section we noted that the geometry 
of some of the clusters with quadratic irrational wind- 
ing exhibit oscillations. It is thus not clear whether they 
have the same fractal dimension D. In this subsection we 
provide numerical test of the claim that the quadratic ir- 
rationals belong to the same universality class. In the 
the following sections we address this question using ad- 
ditional tools. 

To study quantitatively the oscillatory fractal geom- 

In 



ctry we consider the dependence of F[ n ^ on n 
Fig. [impanel a we present compensated plots of Fx (\/2) 

vs.F{^(v / 3) and vs. F± of a DLA as a function of n. 

It appears that although this ratio exhibits oscilla- 
tions, these are bounded and decreasing in amplitude, 
at least up to n — 10 5 . For comparison we show in panel 
b of Fig. |l| a plot of F[ n \y/2)/F[ n) (DLA). Here we 
see the clear difference in dimension as seen in the ra- 
tio approaching zero as a power law in n. In Fig. 

(n) 

we show compensated plots of i<\ of the clusters in 
Figs. |7- 11 versus F± of the golden mean growth. We 
see oscillations on the logarithmic scale, but again these 
are bounded, and we propose that this points towards the 
possibility that all quadratic irrationals winding numbers 
lead to the same overall dimension of the cluster. In the 
next section we address the issue of universality classes 
using additional tools. 



V. TOWARDS A SCALING THEORY: WINDING 
WITH RATIONAL APPROXIMANTS 

To gain understanding of the geometry of the clusters 
grown with quadratic irrational winding numbers we will 
make use now of the well known fact that these irrationals 
can be systematically approximated by rational approxi- 
mants. Thus, having a cluster constructed with a golden 
mean itinerary, a natural question is what happens to the 
growth pattern when p is replaced by ratios of successive 
Fibonacci numbers which are defined by the recursion 
relation F m+1 = F m + F = Q, Fx = 1. Using 



rational approximants p m — F m -i/ F m , the itinerary be- 
comes periodic on the unit circle with period F m , and 
it is observed in simulations (see Fig. || ) that while for 
small clusters n <C n c (F m ) the cluster appears fractal, 
for n ^> n c (F m ) the cluster consists of a set of F m rays, 
sometimes fused into a smaller set of 1-dimensional rays 
whose number is extremely sensitive to the initial condi- 
tions (here controlled by the value of Ao). 



A. The 1-dimensional phase 

The properties of the 1-dimensional phase are impor- 
tant for developing a scaling theory. As an example of the 
interesting behaviour seen as a function of Ao consider 
Fig. [l^ in which clusters with W — 144/89 are grown 
with 4 values of A which are 0.11, 0.22, 0.44, and 0.88. 
Evidently the cross over from fractal to 1-dimensional 
behaviour depends on Ao- We also note that the number 
of rays in the 1-dimensional phase has a nonmonotonic 
dependence on Ao. This indicates high sensitivity of the 
number of rays to changes in the initial conditions. Obvi- 
ously, the radius of the cluster in the 1-dimensional case 
is inversely proportional to the number of rays. On the 
other hand, we have found a surprising invariant: F^ 
is asymptotically invariant to the number of rays (i.e. to 
initial conditions) being always equal to n\/~\o/Q, up to a 
constant of proportionality depending on the microscopic 
parameter a only. The numerical evidence is shown in 
Fig. [n]. Note the convergence to the golden mean in 
panel a, and to \[2 in panel b (which is the value of 
the ratios of \/Ao)- This finding puts strict bounds on 
the number of possible rays. The upper bound is ob- 
viously Q. The lower bound stems from the inequality 
Rn < AF^ , meaning that the number of rays must be 
larger than Q/4. This invariance also indicates that the 
geometry of the rays is not arbitrary, and that the angles 
between them are arranged to agree with an invariant 



(n) 



B. Scaling Function 



The crossover in fractal shape is a general result for 
any periodic itinerary with W = P/Q, and suggests the 
existence of a scaling for F± of the form 



(18) 



where we have assumed that the crossover cluster size 
scales as 



n c (Q) ~ Q c 



(19) 



The asymptotic forms of F{ obey F^ ' ~ n 

for n <C n c (Q), while F[ n ^ ~ {n/Q)^/X{) for n > n c (Q). 
In the first asymptote we expect D to be the same for 
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all values of rational approximants to p, including the 
limiting fractal cluster. The growing cluster cannot dis- 
tinguish between the rational approximant and the lim- 
iting irrational as long as the fractal phase is observed. 
The second asymptote is demonstrated in the previous 
subsection. Thus we require that the asymptotic forms 
of the scaling function obey 



f(u) -> /(0) as u -> , 
f(u) ~ u as u — > oo . 



(20) 
(21) 



The second asymptote (|2l]) determines the scaling rela- 
tion 



(22) 

1.86 and consequently 



a = D/(D - 1) 

For the Golden Mean fractal D ss 

in this case a « 2.16. 

? (n) 



In Fig. HI F[ n) /{n 1 ' ^) is plotted against the scal- 
ing variable u = n x / a /Q for six different clusters with 
different values of W and Aq. The best data collapse 
was obtained using the value a = 2.15 The data col- 
lapse achieved is readily apparent with the scaling func- 
tion f(u) predicted by the theory. 



VI. THE CROSSOVER AND THE ESTIMATE OF 
THE DIMENSION 

In this section we discuss the properties of the confor- 
mal map 4>\£ which determine the cross over from fractal 
to 1-dimensional growth. In other words, we will attempt 
to provide an independent estimate of function 
of the winding number W . If we succeeded to estimate 
the exponent a in (19) independently from Eq.(p2]), we 
would have an equation for the dimension. 

To understand the crossover, we note that the reason 
for the fractal growth phase with rational winding is that 
after every event of growth the interface z(e 10 ) is non- 
locally reparametrized in addition to the local growth 
event. Accordingly, a periodic orbit on the unit cir- 
cle is not necessarily mapped to a periodic orbit in z. 
The region in the unit circle which is significantly af- 
fected by growing the nth bump has a scale \f\^, cen- 
tered around n |l5| . Accordingly we can estimate when 
reparametrization will cause a "miss" in the mapped or- 
bit: as long as 



An > 



2tt 

~Q 



(23) 



the growth will remain fractal. We can therefore expect a 
crossover to 1-dimensional growth when this condition is 
violated, something that is bound to happen since typical 
values of A„ are expected to decreases with n, cf. Eq.jT^) 
and the discussion below. 

What remains is to estimate A„ as a function of n in 
the cross over region that is defined by 



2n/Q 



(24) 



In the fractal region A„ is a highly erratic function of n. 
Even though we do not have here randomness in the sense 
of DLA, it is natural to consider, in a fashion similar to 
Eq. ( |l4[ ) , the distribution of A^ over Q successive steps of 
growth. For Q large enough such distributions have well 
defined moments. In particular consider the first moment 



1 - 



(25) 



k—n — Q 



The power law dependence of and Eq.(|To|) imply 

that this moment has to be 



(A„)q 



anD 



(26) 



If we estimate A„ c in Eq. ( |24| ) by its mean (pq), we 
would write 



A„ c ~ l/n c 



nc 



Q 2 . 



(27) 



Thus D/(D - 1) = 2 or D = 2. Even though we get an 
overestimate, this is a good indication that we are on the 
right track. The reason for the overestimate is that we 
neglected the fluctuations that sometime lead to A n much 
larger than the mean. We expect a cross over to occur 
when the largest VAfc are smaller than 2ir/Q 2 , since it is 
enough to have a few large Afc to cause a reparametriza- 
tion that will ruin a potential periodic orbit. We thus 
seek a condition 



A- ax = max{A fe }]U-Q « ^ • 



(28) 



We note that A^ is an erratic function of fc, and therefore 
the condition d2S|) ca n be met more than once in a given 
series A&. In Fig. fill we show two log- log plots of n c com- 
puted from the value of n for which */ A" lax = 2n/Q, plot- 
ted as a function of Q = F m . The cross over value n c was 
computed in two different ways. In circles we exhibit the 
values obtained from measuring when ^/A™ 3 * = 2n/Q 
for the first time, and in squares we exhibit the values 
obtained from yj A™ ax = 2n/Q for the last time. Com- 
puting the slopes by linear regression and averaging be- 
tween them we find the scaling law 



Q 



2.17±0.03 



(29) 



Comparing with Eqs.([T^), ( p2|) we get an estimate for 
D = 1.86 ± 0.03, in excellent agreement with the deter- 
mination of the dimension by . 

We note in passing that A™ ax can be assigned a gen- 
eralized dimension in the language of Hentschel and 
Procaccia Eql. Define 



i 

Q 



E 



(30) 
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From pCI the precise scaling law is 



lim (A 



1) 1/q ~ n- 2D °°/ D 
n/Q n 



(31) 



Comparing with ( p2[ ) we conclude that in this case there 
exists a scaling relation 



D - 1 



(32) 



Such a scaling relation was conjectured by Turkevich and 
Scher for DLA 17 (of course with a different D and Doc). 
While there are severe doubts about the correctness of 
this conjecture for DLA |Q, we point out that in our 
case it follows directly from elementary considerations. 



A. The period doubling itinerary 

Even though the period doubling itinerary leads to a 
cluster whose fractal dimension differs from the quadratic 
irrational windings, we show here that the ideas pre- 
sented above pertain equally to this growth pattern. In- 
stead of rational approximants we use here, naturally, 
2 n -periodic orbits which are obtained by cutting the 
itinerary ( |l7| ) after 2™ iterations and repeating it periodi- 
cally. The crossover from fractal to 1-dimensional growth 
is seen also in this case, and we can use it in a very similar 
way to identify the crucial exponent that determines the 
dimension of the asymptotic cluster. Indeed, the whole 
set of ideas developed above repeats verbatim by chang- 
ing Q with 2™ . What remains is to find A™ ax as a function 
of n. In Fig .|2^ we show the data collapse obtained as 
in Fig. |l8| for the quasiperiodic analog. We show nine 
different data sets with periodic itineraries of periods 32, 
64 and 128 and y% values of 0.22, 0.44 and 0.88. The 
scaling function for these data sets is plotted as a func- 
tion of u — n A4: /Q, where the exponent is computed from 
D = 1.78. It is noteworthy that the scaling function ob- 
tained appears identical to the scaling function f(u) for 
the quasiperiodic family. For comparison we added in 
Fig .go] also one curve from the quasiperiodic class, and 
it appears indistinguishable from the rest. 

The conclusion from this data collapse is that the 
mechanism governing the crossover from fractal to 1- 
dimensional growth phases here is the same as the one 
discussed above for the quasiperiodic itineraries. The 
difference between the dimensions of the period doubling 
cluster and the quasiperiodic cluster must lie in the differ- 
ent numerical value of the exponent characterizing A™ ax 
as a function of n. In this case the natural averaging cy- 
cles are of length Q = 2". Fig. |2l| is the analog of Fig. |l| 
for the period doubling itinerary, where the critical value 
n c was estimated from the first time that U A™ ax became 
smaller than 2ir/2 n . 



Computing D we find D = 1.75±0.05 in good agreement 
with the numerical estimate from . 



VII. SUMMARY AND THE ROAD AHEAD 

The main points of this paper are as follows: 

• The iterated conformal maps algorithm for fractal 
growth patterns offers a convenient way to intro- 
duce a large number of deterministic growth models 
with highly non-trivial fractal geometry. 

• Itineraries with irrational winding numbers gener- 
ate fractal growth patterns. We proposed that all 
the quadratic irrationals produce clusters of the 
same fractal dimension, in spite of different appear- 
ance. 

• By considering a series of rational approximants 
we could produce a scaling theory of the growing 
clusters, achieving data collapse for all values of 
n, Xq and P/Q- 

• Identifying the mechanism for the cross over from 
fractal to 1-dimensional growth phases we could 
pinpoint the exponent that determines the fractal 
dimension D. This exponent characterizes the n 
dependence of the extremal values of X n . 

• The mechanism appears general; itineraries leading 
to different cluster dimensions, like the period dou- 
bling itinerary (|l7]) and its truncated versions, can 
be understood in the same way. The scaling func- 
tion (|l^) and the scaling relation (|2^) are general, 
but the exponent a changes. Its determination by 
the scaling of A™ ax Eq.(|2^) is however general. 

We note that all the numerical tests point out in favour of 
this scenario, and in our opinion rule out a value D = 2 
for the clusters discussed above. The only way to get 
2-dimensional growth, as shown above, is if the distribu- 
tion of X n does not multiscale, i.e. all D q are the same, 
and the scaling of A™ ax identifies with the scaling of the 
average of A n . 

Nevertheless, we point out that the crucial step in our 
scenario, the determination of the exponent a in Eq.(fl9|), 
was achieved numerically. The scaling theory presented 
above has a strong flavour of a renormalization group ap- 
proach. It appears that such an underlying theory may 
have a low codimension, maybe with 1 important expo- 
nent, the one characterizing the rate of crossover of the 
rational approximants to the irrational limit. The search 
of such a theory appears to be an important task for the 
near future. 



The linear regression provides us with the the scaling 
law 



i2.33±0.1 



(33) 



G 



ACKNOWLEDGMENTS 

Wc benefitted from discussions with T.C. Halsey, 
Trcsser and L. Peliti. This work has been suppor 
in part by the European Commission under the T] 
program and the Naftali and Anna Backcnroth-Bron 
Fund for Research in Chaos and Complexity. 



[1] TA. Witten and L.M. Sander, Phys. Rev. Lett, 47, 1 
(1981). 

[2] L. Niemeyer, L. Pietronero and H.J. Wiessmann, P] 

Rev. Lett. 52, 1033 (1984). 
[3] L. Paterson, Phys. Rev. Lett. 52, 1621 (1984). 
[4] R. M. Brady and R. C. Ball, Nature 309, 225 (1984); D. 

Grier, E. Ben- Jacob, R. Clarke, and L. M. Sander, Phys. 

Rev. Lett. 56, 1264 (1986) 
[5] A. Erzan, L. Pietronero, A. Vespignani, Rev. Mod. PI 

67, 545 (1995). 
[6] P. Meakin, in Phase Transitions and Critical Phenomi 

eds. C. Domb and J. L. Lebowitz, Vol.12 (Academic, ] 

1988). 

[7] P.G. Saffman and G.I. Taylor, Proc. Roy. Soc. Lon 

Series A, 245,312 (1958). 
[8] B. Shraiman and D. Bensimon, Phys. Rev. A30, 2 

(1984) . 

[9] M. B. Hastings and L. S. Levitov, Physica D 116, 
(1998). 

[10] B. Davidovich, H.G.E. Hentschel, Z. Olami, I.Procac 
L. M. Sander, and E. Somfai, Phys. Rev. E, 59 1 
(1999) 

[11] B. Davidovich and I. Procaccia, Euro. Phys. Lett, 
547 (1999) 

[12] T.C. Halsey, P. Meakin and I. Procaccia, Phys. Rev. L 
56, 854 (1986). 

[13] P. L. Duren, Univalent Functions, (Springer- Verlag, New 
York, 1983). 

[14] T.C. Halsey, Phys. Rev. Lett. 59, 2067 (1987). 
[15] T.C. Halsey, Phys. Rev. Lett. 72 1228 (1994). 
[16] H.G.E. Hentschel and I. Procaccia, Physica D 8, 435 
(1983). 

[17] LA. Turkevich and H. Scher, Phys. Rev. Lett. 55, 1026 

(1985) . 




FIG. 1. A DLA cluster, n = 10 5 . 




FIG. 2. A Bernoulli cluster, n = 10 5 . 
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FIG. 4. The average of kXk over the last M iterations 
of the Bernoulli itinerary, with M — 10 3 (continuous line) 
and M = 10 4 (dashed line). The horizontal line is the ex- 
pected value, 1/aD. The fluctuations are typical, reflecting 
the multi-scaling distributions of in which large deviations 
are highly probable. 




FIG. 6. The cluster grown with W = p to n = I0 5 
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FIG. 12. The cluster grown with the period doubling 
itinerary to n = 10 5 FIG. 14. Panel a: The ratio between of F^ n) of the clusters 



grown with W = 2 1/2 and W = 3 1/2 . Panel b: The ratio 
between of F[ n ^ of the cluster grown with W = 2 1/2 and a 
typical DLA 
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FIG. 15. The ratio between of i*\ of the golden mean 
and other quadratic irrationals. Panels a-d show respectively 
W = 2 1/2 , W = (1 + (10) 1/2 )/3 , I¥((13 1/2 - l)/2) and 
W = 7 1/2 



FIG. 16. Clusters grown with W = 144/89 with four dif- 
ferent values of Ao, from 0.11 to 0.88. 
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FIG. 17. Panel a: The ratio of F{ n} for clusters grown with 
A = 0.88, and winding numbers W = 89/55 and W = 144/89 
(upper curve) and W = 144/89 and W = 233/144) (lower 
curve). Both converge to the golden mean. Panel b: Sim- 
ilar plots with W — 144/89. Ao = 0.22 is compensated by 
Ao = 0.11, 0.44 by 0.22 and 0.88 by 0.44. All these plots 
converge to 2 1 / 2 . 
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FIG. 18. Scaling behaviour for six separate data sets for 
f(u) with u = n A6 /Q. Shown are AO = 0.88 with W = 89/55, 
W = 144/89, W = 233/144, W = 377/233, AO = 0.44 with 
W = 144/89, AO = 0.22 with W = 144/89. 
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FIG. 20. Scaling behaviour for nine separate data sets for 
clusters grown with truncated period doubling itineraries, in 
addition to one data set of the quasiperiodic class. See text 
for details 
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FIG. 19. The measured the cross over values n c as a func- 
tion of Q in a log-log plots. In dots are the values of n for 
which the condition A™ ax = Atv 2 /Q 2 was met for the first 
time, in squares for the last time. 
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FIG. 21. The measured the cross over values func- 
tion of Q — 2™ in a log-log plot. Shown are the values of n 
for which the condition A™ ax = (27r/2 n ) 2 was met for the first 
time. 
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